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ABSTRACT 

Extremely luminous, red eruptive variables like RV in M31 are being suggested as 
exemplars of a new class of astrophysical object. Our greatly extended series of nova 
simulations shows that classical nova models can produce very red, luminous eruptions. 
In a poorly studied corner of 3-D nova parameter space (very cold, low-mass white 
dwarfs, accreting at very low rates) we find bona fide classical novae that are very lu- 
minous and red because they eject very slowly moving, massive envelopes. A crucial 
prediction of these nova models - in contrast to the predictions of merging star ("merge- 
burst") models - is that a hot remnant, the underlying white dwarf, will emerge after 
the massive ejected envelope has expanded enough to become optically thin. This blue 
remnant must fade on a timescale of decades - much faster than a "mergeburst" , which 
must fade on timescales of millenia or longer. Furthermore, the cooling nova white 
dwarf and its expanding ejecta must become redder in the years after eruption, while a 
contracting mergeburst must become hotter and bluer. We predict that red novae will 
always brighten to L ~ IOOOLq for about 1 year before rising to maximum luminosity 
at L ~ 10 6 - 1O 7 L . The maximum luminosity attainable by a nova is likely to be 
L ~ 1O 7 L0, corresponding to M ~ —12. In an accompanying paper we describe a 
fading, luminous blue candidate for the remnant of M31-RV; it is observed with HST 
to be compatible only with the nova model. 

Subject headings: accretion, accretion disks — binaries: close — novae, cataclysmic 
variables — white dwarfs 



1. INTRODUCTION 

The realization that classical novae are colossal explosions of previously nondescript stars 
was hammered home to the astronomical community with the appearance of first-magnitude GK 
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Persei — Nova 1901 A.D. Outshining all but a handful of stars in the night sky, it jump-started 
the scie ntific study of classical nova eruptions (jPavne-Gaposhkin Ill957l). The ex traordinary light 
echoes ( Ritchev 1901 ; Couderc 1939 ). and ejected shells of matter (|Ritchev~ 1917 ) associated with 
GK Per are the obvious hallmarks of a highly energet ic, mass-ejecting event. The detection of 
erupting c lassical novae in the Small Magellanic Cloud ( McKibbenlll95ll ). Large Magellanic Cloud 
(|Gilllll927l ) and in M31 (lHubblejll929l ) established novae as ubiquitous beacons radiating at tens of 
thousands of times the luminosity of the Sun or more — a new class of astrophysical phenomenon 
( Harwiti 1 19811 ). Over the past century, of the order of a th ousand erupting class ic al novae have 



been detected in the Galaxy, in gala xies of the Local Group (ICiardullo et al.lll990l ). (jShara et al 



2004 ) and as far awa y as the Virgo (jFerrarese et al.ll2003l : IShara et al. 1120041 ) and Fornax clusters 
dNeill Sharall2004h . 



Less appreciated but just as important is that these early nova detections demonstrated the 
existence of close, interacting binary stars in galaxies outsid e our own. This could be appreciated 
only after two more seminal discoveries. First. I Walker I (|1954j) demonstrat ed that nova Herculis 1934 
(DQ Her) was a short period (4.65 hour) binary. A decade later, Kraft's (|Kraftlll964l ) spectroscopy 
demonstrated that old novae must be composed of a white dwarf (WD) accreting hydrogen-rich 
matter via an accretion disk from a red dwarf (RD) companion. Two logical consequences of this 
model are that (1) classical nova e are powered by therm onuclear runaways (TNR) in their white 
dwarfs' hydrogen-rich envelopes (|Starrfield et al. Ill972i) and (2) novae must self-extin guish when 
their erupting envelopes are all but completely ejected (jPrialnik. Shara k, Shavivlll977l ). 



Shortly afterwards it was realized that three parameters (WD mass, WD temperature or 
lumin osity, and WD mass accretion ra t e) determine t he physical characteristics of all nova erup- 
tions (jShara. Prialnik. &: Shaviv I Il980l ; iPrialnikl Il995l ) . Increasingly sophisticated input physics 
and rapidly increasing computer power led to self-consistent nova models including accretion, ther- 
monuclear runaways (TNR) and ejection o yer many nova cycl e s. Th ese culminated in the extensive 
grid of multi-cycle nova eruption models (jPrialnik Koveta (|1995l ). hereafter Paper I). Covering 
what was then thought to be the entire 3-D parameter space that would produce novae, this set of 
models reproduce d the range of energetics, timescales and mass ejections attributed at that time 
to classical novae (IWarner 1 11995). 



So me remarkable red, luminous eru ptive variables are often claimed (e.g. in many of the pa- 



pers m iCorradi. R. &: Munari. U.I (|2007l )) to have physical properties inconsistent with all of the 
nova models in Paper I. While thermonuclear-powered nova eruptions can explain the luminosi- 
ties of these objects , their massive and very co ol ejecta seem t o tally di fferent from those o f othe r 



novae 



dMould et al.l (Il99qh bond et al.1 (hooal ). 



Munari et al 



(1200211. baneriee fc Ashokl (hood) 



Kimeswenger et al. I (|2002h . ISoker fc Tvlenda I (120031 ) JKipper. Klochkova. h Annuk I (120041 ) . iLvnch et al 
(|2004l )). Does the nova model fail for these objects? Does this signal the existence of a hitherto 
unknown stellar eruptive phenomenon? Although V838 Mon appears to be incompatible with a 
nova model (as detailed below), are the other luminous red variables similarly constrained? This 
oft-repeated claim that M31-RV and similar, ultra-luminous red, eruptive variables cannot be ex- 
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plained as classical novae has prompted us to reexamine the 3-D parameter space that produces 
classical novae. The goal of this paper is to determine whether very luminous, red eruptive variables 
can, or cannot be produced by WD thermonuclear events. We also want to determine whether there 
are one or more observational signatures that can clearly differentiate between thermonuclear novae 
and other models of very luminous, eruptive red variables. If such red novae emerge naturally from 
our simulations then we might expect them to appear in nature. Of course, success in these efforts 
will not prove that thermonuclear novae are responsible for M31-RV or any given similar objects. 
But it will demonstrate that thermonuclear nova models should not be dismissed out of hand as the 
explanation for a very luminous red variable unless other evidence excludes such an interpretation. 
Only in the case of V838 Mon does such evidence seem to be in hand. 

A pioneering attempt to consider very col d, low mass whit e dwa rfs accreting hydrogen to 



produce very massive ejected shells was made by llben h Tutukovl ()1992l ). Such massive ejecta are 



essential to explain the cool red spectra. We have extended our grid of multi-cycle novae to lower 
white dwarf masses, lower accretion rates and colder white dwarfs than ever before, seeking the hard 
boundaries of nova 3-D space. We have found many new combinations that do give rise to classical 
novae, including those with t he most ext r eme l uminosities and ejected masses yet encountered. 



These models are described in lYaron et al.l (|2005l ) (hereafter Paper II). Here we concentrate on two 



of those models and several additional ones that yield better fits to the basic characteristics of M31 
RV than has been achieved before by any nova simulation. In an accompanying paper we report 
the discovery of a luminous, very blue star very close to the location of M31-RV that is behaving 
like the remnant of a classical nova eruption. Over the past decade that star has become much 
redder and fainter, and it now resembles old classical novae. 

In § [2] we summarize the key observations that guide our new nova models. In § [3] we detail 
the initial models, and in § H] we describe their evolution during nova eruptions. We compare the 
observed characteristics (detailed in § [2]) with the models in § [5j and offer testable predictions con- 
cerning M31-RV and the M85 eruptive variable. We briefly summarize our results and conclusions 
in §E3 



2. OBSERVATIONAL CONSTRAINTS 



2.1. M31 RV 



Determining accurate distances, and hence luminosities, for Milky Way stars is observationally 
challenging. It is no accident that classical novae and Cepheid variables were first calibrated as 
standard candles in the Magellanic Clouds and in M31. The well-established distance of the An- 
dromeda galaxy makes the relatively well-studied system M31- RV ("RV" stands for red variable) 
the defining extragalactic member of its class. As pointed out bv lRich et al.l (|1989l ). at peak bright- 
ness M31-RV approached a bolometric magnitude Mb i = —9.6, making it the most luminous red 
supergiant in the Local Group. The key observational constraints on this remarkable star have been 
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presen ted b ylRich et al.l ( 1989 



Mould et all |l990|), Bryan fc Roverl (11992 ). iTomanev fc Shafter 



(|1992l ). and lBoschi fc Munaril (|20Q4| ) . These observed characteristics and constraints (summarized 



below) must be mimicked by any model that seeks to explain this and similar outbursts. 



6. 



The maxim um observed brig htnesses attained by M31 RV were V=15.3 (jBryan &: Royerlll992l ) 
and 1=14.9 ( Rich et al.lll989 ). The maximum absolute magnitude was thus at least M = —9.3, 
or 4 x 1O 5 L0. The emission-line spectrum, together with the unknown dust absorption, 
precludes a definitive measure of the absolute magnitude at maximum. Using reasonable 
assumptions abo ut reddening, M g iant bolometric corrections and an adopted M31 distance 
modulus of 24.2. iRich etall (|l989h deduced M bol = -9.96, corresponding to 8 x 10 5 L Q . The 
other authors noted above give values similar to this for the peak outburst luminosity of M31 
RV. 



2. The progenitor of M31 RV was fainter than I = 20.4 (jMould et al.lll990l ). thereby demon- 
strating an outburst amplitude exceeding 5.5 magnitudes. 

3. The time to brighten the last 3 or 4 magnitudes t o maximum light was less than 29 days, and 
could have been as short as hours or a few days (jBryan Roverl ll992). 



M31 RV faded to 6 x 10 4 L Q at an epoch 70 days after maximum observed brightness. It 
faded below 10 4 Lq 290 days after maximum brightness. 

M31 RV displayed very red colors in the Gunn system (g-r = 1.33) and in the Johnson system 
(B-V = 1.9 and V-R = 1.0) near maximum light. At the same time it displayed a spectrum 
resembling that of an MO supergiant. Infrared observations in the 1 to 3.6 micron range 
demonstrated that the ejected shell became as cool as 4000 Kelvins at maximum light, with 
a photospheric radius of 2000-R©. About 70 days later the effective temperature dropped to 
1000 Kelvins with a photospheric radius of 8000-R©. 

The "coronal" or nebular phase of emission, seen in virtually all classical novae, did not appear 
in M31 RV. An ejected shell mass considerably larger than that of most novae therefore seems 
inescapable. 



7. The veloci t y of e jection of material from M31 RV was ~ 150 km s 1 ( Mould et al.lll990i ) 



Rich et al.l (|1989l ) reported a +390 km s 1 radial velocity c hange in the variable over a span 



of eight hours, in both Ha and stellar absorption lines. ITomanev Shafter I (|1992l ) saw 
no such change, but did report significant [Nil] emission lines flanking Ha, which might be 
rapidly variable. 
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2.2. The M85 Optical Transient 

This object is the most luminous of all the "red novae". 



1. At a distan ce of 15 Mpc the ab solute R-band magnitude of this object at maximum was 
close to -12 (Kulkarni et al.1 120071 ). At 180 days after maximum it was producing 3 x 1O 5 L0 
(|Rau et alll2007h . 



This "red nova" brightened by at least 7.7 magnitudes from its pre-eruption state. 

The object faded most rapidly in blue light, and most slowly in the red. Nearly peak brightness 
was attained for 70 days, and the fading time thereafter was a few months. 



4. Optical spectra show that the shell effective temperature was 4600 Kelvins near peak lu- 
minosity. Infrared multicolor images determine the ejecta's effective temperature to be 950 
Kelvins at 6 months after the eruption began. 

5. The optical material ejected from the M85 optical transient displayed a FWHM velocity of 
350 km s , and an infrared expansion velocity of 870 km s _1 . 

6. Like M31-RV, t he M85 optical tr ansient (M85- OT) was not associated with any bright star 
forming region (jOfek et al II2008I ). The g- and z-band absolute magnitudes of the progenitor 
were fainter than about -4 and -6 mag, respectively, corresponding to an upper limit for a 
progenitor (main sequence) mass of 7Mq. 

After this paper was co mpleted a similar lum inous red eruptive variable was discovered in the 
the Virgo galaxy M99 (jKasliwal et all 12010 ) . Hubble Space Telescope archival imagery rules 
out red supergiants, blue supergiants and single main sequence stars earlier than type B2 as 
the progenitor. 



2.3. V4332 Sgr 



The eruptive va riable V4332 Sgr be ars remarkable similarities to the luminous red objects we 
have already noted ([Martini et allll999l ). There is not, unfortunately, a well-determined distance 
to V4332 Sgr, and hence its luminosity is quite uncertain. However, it is the only object of this 
class to have photometric observations in the years just b efore its eruption. A m inimum fourfold 
brightening of V4332 Sgr in the decade before its eruption (jGoranskii et al. 1120071 ) is an important 
constraint on models of this star. 
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2.4. V838 Mon 



V838 Mon riveted the attention of astronomers when it erupted in 2002 January for two 
reasons. First, it nearly reached naked-eye brightness, peaking at V = 6.7 . Even more remarkable 
were the spectacular series of light echos surroundi ng the object. H igh resolution Hubble Space 
Telescope imagery and polarimetry of those echoes (|Bond et al.ll2003l ) conclusively established the 
distance of V838 Mon at 6 kpc or greater. This in turn has allowed strong constraints on the 
physical characteristics of the variable, particularly its luminosity. At a distance of at least 6 
kpc (and reddenin g E(B-V) = 0.8) the absolute magnitude of V838 Mon at maximum was at 



least M v = —9.6 (jBond et al.1 120031 ). V838 Mon and M31 RV are thus remarkably similar in 
their overall energetics, as well as their spectroscopic and photometric behaviors. The excellent 
distances available for each object have been crucial in demonstrating this remarkable agreement. 
The extremely red colors and cool spectra in the decline phase demand ejecta much more massive 
than has hitherto been predicted by nova models. 

There is strong evidence t hat V838 Mon is asso ciated with a B3V dwarf star, which is itself a 
member of a young association (lAfsar &; Bondll2007l ) of at least 4 B-type stars. Such an association 
is too young to produce a cataclysmic binary, especia lly with a low mass white dwarf, unless 
strong dynamical stellar interactions and collisions occur (jPortegies Zwart et al. Il2010i ) . (A possible 
formation mechanism would have the CO white dwarf core of a relatively massive star liberated 
from its envelope via collision with a neutron star, with the white dwarf capturing the neutron 
star's low mass companion in the process. But even then the time to cool the white dwarf, and to 
slowly accrete a massive envelope is much longer than the lifetime of a B3V star) . Ordinary binary 
star evolution also cannot form a nova progenitor for V838 Mon within the main sequence lifetime 
of the B3V star. 

Ther e appear not to b e correspondingly young stellar associations at the sites of M31-RV or 
M85-OT (jOfek et al 1 120081 ) or M99-OT, so that very luminous, red classical novae like those we 
describe below remain viable models for these events. 



3. NEW NOVA MODELS 

In Paper II we published the results of evolutionary calculations of nova outbursts for a wide 
range of parameter combinations, spanning the 3-dimensional parameter space of novae: white 
dwarf mass Mwd, white dwarf core temperature Ty/D (correlated with the WD luminosity and 
hence age) and accretion rate M. This grid of models was an extension to the grid published in 
Paper I ten years earlier. The parameter values adopted in the extended grid were: four My/D 
values— 0.65, 1.00, 1.25 & 1.40 M , three T WD values— 10, 30 & 50 x 10 6 K and eight M values— 
10~ 6 through 10~ 12 M yr" 1 , as well as 5 x 10~ 13 M yr -1 . 

The hydrodynamic Lagrangian stellar evolution code used in all our studies is described in 
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some detail in Paper I. It includes OPAL opacities (jlglesias fc Rogers 1119961 ). an extended nuclear 
reactions network comprised of 40 heavy element isotopes, and a mass-loss algorithm that applies 
a steady, optically thick supersonic wind solution (following the phase of rapid expansion). In 
addition, diffusion is computed for all elements, accretional heating is taken into account and con- 
vective fluxes are calculated according to the mixing length theory. Initial models were prepared 
for the four WD-mass values and three temperatures by cooling WD models from higher temper- 
atures. Each nova model was followed through several consecutive outburst cycles. In each case, 
one cycle was then chosen as representative. In Paper II we showed — by analytical considerations 
supported by numerical calculation results — that the parameter space where nova outbursts oc- 
cur is limited. In order to test the extremes, we added calculations for still lower accretion rates 
(5 x 1CT 13 Mq yr -1 ), a few cases for Mwd = 0AM & , and several cases for T\yd of only a few 
x 10 6 K. Many of these cases were found to lie outside the classical nova outburst parameter space; 
no eruptions ever occurred. Others presented features that appeared uncharacteristic of observed 
classical novae, but rather reminiscent of the properties of M31 RV. The results of two such extreme 
cases of the Paper II grid are given in Table [H where m acc is the accreted mass, m e j - the ejected 
mass, Z e j - the mass fraction of heavy elements in the ejecta, v av - the average expansion velocity, 
Tmax — the maximal temperature attained in the burning shell at the base of the envelope, L max 
- the peak luminosity, and t^^oi ~ the time of decline of the bolometric luminosity by 3 mag. In 
the case of these very cold WDS, consecutive cycles are not identical, while they were identical for 
the hotter WDs considered. The differences between cycles are due to heating of the WDs by nova 
eruptions (see below). The results shown are for the first cycle, not as representative of all such 
eruptions, but as a reasonable example. We emphasize these models' unusual combination of very 
large ejected masses, low ejection velocities and very high (super-Eddington) peak luminosities. 

M31 RV must have ejected a shell in excess of 1O" 3 M to have become as cool as was observed. 
A very cold low-mass WD, accreting at a low rate is clearly the candi d ate to look for, if a very high 
ejected mass is sought, as was also pointed out by llben Tutukovl (|1992l ). However, the O.4M 



WD eruption declines too slowly to match the observations of M31-RV. We therefore decided to 
adopt an initial mass of 0.5M Q , which is also more realistic for a C-0 WD. We allowed the WD 
model to cool to core temperatures of 4 x 10 6 K (3.4 x 10 9 yr), 3 x 10 6 K (5.1 x 10 9 yr) and finally, 
2 x 10 6 K (8.8 x 10 9 yr). For each model we found the lowest accretion rate that still produced nova 
outbursts. For the middle T\yd temperature, we also ran a model with a higher accretion rate, for 
comparison. Each model was evolved through several nova outburst cycles, as in Papers I and II. 
The results of these calculations are given in Table 

In contrast to the previous models, for which the WDs were hotter, the present, very cold 
models have very long cooling timescales, and thus the outer layers of the WD, which are heated 
during an outburst, do not cool back to the initial temperature before the next eruption. This is 
illustrated in Figure 1, where the new models are compared with a series of models from Paper II. 
Consequently, consecutive outbursts are not identical, but rather change monotonically, although 
slowly. Since the accreted mass decreases with increasing temperature of the outer layers, we list 
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in Table [2] the results for the first outburst in each run, which produces the highest ejected mass. 

In Table [3] we list abundances (by mass) of the main CNO isotopes, as well as isotopic ratios. 
We draw attention to the high oxygen content of the ejecta relatiye to carbon. This is typical of 



nova outbursts on low- mass WDs, as shown by lKovetz fc Prialnild (119971 ): the trend is reversed for 
massive WDs. 

Some of the parameters' trends reported in Tables 2 and 3 are not mo notonic; the explan ation 



is as follows. The nova phenomenon is a 3-parameter family of events (jYaron et al.l 120051 ) . and 
dependence on parameter values is not linear. For a given WD mass and accretion rate, Z decreases 
with decreasing WD temperature. This is because more mass is accreted before the outburst when 
the WD is cold, while diffusion of hydrogen into the core at the inner boundary of the accreted 
shell slows down. As a result, the mixture of core material into the envelope (once convection sets 
in at ignition) is more diluted. This explains the differences between the first two rows in Table 2. 

For given WD mass and temperature, Z increases as M decreases because the accretion time 
is considerably longer and there is more time for diffusion. However, eventually, diffusion reaches 
equilibrium and effectively ceases and this is independent of M. Therefore, if M decreases further, 
Z will decrease due to dilution. This explains the difference between rows 4 and 2 in Table 2. When 
two of the three parameters are changed at once, and these parameters each tend to change Z in 
the opposite sense, it is even more difficult to predict the outcome. This explains why the results 
of line 3 in Table 2 may not be interpolated from the other results of the table. 

In order to test the lowest possible core temperatures and their effect on the nova charac- 
teristics, we cooled three low-mass WD models for approximately a Hubble time (1.3 x 10 10 yr) 
and then started accretion. At such low temperature an accretion rate of 1O -11 M0 yr -1 did not 
produce a thermonuclear runaway (see Paper II for a discussion of the limited parameter space for 
TNR). We therefore adopted M = 1O -1O M0 yr -1 for all these models. The results — again, for 
the first cycle of the series — are summarized in Table 4. 

Generally, the results are similar — high ejected masses, low velocities and very high peak 
luminosities — although they differ in detail. Nova outburst characteristics a re sensitive to even 



modes t variations of one of the three basic parameters (see Papers I and II, and ISchwartzman et al 



()1994l )). Since small effects may arise from factors that are not included in the model (e.g. rotation 
or the presence of the secondary star) we do not aim for a precise match of any model with a 
particular ultra-luminous nova. That we are able to achieve good matches to the observations 
without any "fine-tuning" speaks for the robust nature of the models. 

The same arguments about non-motonically behaving parameters apply to Table 4. The mass 
accreted before outburst decreases with increasing WD mass. For a given accretion rate, one would 
expect higher Z with increasing WD mass because the dilution is weaker, but also lower because 
the time available for diffusion is shorter. Moreover, the larger the WD mass, the higher is the 
gravitational acceleration. This makes the diffusion of hydrogen inwards more difficult. It is not 
surprising, therefore, that the change of Z with WD mass is not monotonic. Similar arguments 
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apply to the ejection velocity. Larger WD masses imply higher degeneracy and stronger outbursts 
(with higher maximum temperature), which tend to lead to higher ejection velocities. The higher 
g, however, has the opposite effect on velocity, leading to the non-monotonic behavior of average 
ejection velocity with respect to the WD mass. 

The strong outbursts and high luminosities of these very cold models are the direct result of 
the very high degeneracy attained at the base of the envelope just before the onset of the TNR. As 
an illustration, the Fermi parameter Ef oc (InP — 2.5 In T) ranges from 3.2 to 6.1 for the models in 
Table 2, whereas for the cold models ^Twd = 10 7 K) of Paper II, the Fermi parameter is less than 
2. 



4. MULTIPLE OUTBURST PEAKS: MULTIPLE ERUPTIONS OF 
ULTRA-LUMINOUS RED NOVAE 

The detailed evolution of the coldest O.5M WD model, accreting at a rate of 7x 10 _11 M Q yr _1 
(last but one entry of Table [2]) , is shown in Figures 2 through 7. 

Intense nuclear burning (in highly degenerate matter) occurs for several years before an out- 
burst, as shown in Figures 2, 3 and 4. A year before any detectable optical brightening, a strong 
thermonuclear pulse occurs at the base of the highly degenerate hydrogen-rich envelope, rising to 
10 10 Lq, and lasting only one day. This pulse releases enough energy to expand the white dwarf 
envelope a hundredfold. No subsequent pulse is as powerful, and we emphasize that this pulse, a 
full year before mass-loss commences, provides almost all of the energy needed to unbind the white 
dwarf envelope. 

A very important characteristic of all these models is that mass ejection is not continuous, but 
occurs in a few pulses, separated by several contractions of the remaining envelo pe and subsequent 



rebou nd. This hydrodynamic- nu clear phenomenon was already pointed out by lKovetz &: Prialnik 



(1 19941 ). and discussed in detail by lPrialnik &; Livid (| 19951 ). During the first mass-ejection producing 



pulse, illustrated in Figures 5 and 6, the mass of the nova envelope is sufficiently reduced that 
pressure at the envelope base drops significantly. This, in turn, causes the inner part of the 
envelope to start contracting, and pressure rises at the envelope base again. The enhanced pressure 
obtained in the burning shell at the base of the envelope produces a sufficient increase in nuclear 
luminosity to drive re-expansion of the inner envelope layers. The effect propagates to the surface, 
where a second pulse of mass loss occurs, sometimes characterized by higher expansion velocities. 
In some cases, a third, similar pulse occurs before the remnant envelope mass becomes sufficiently 
low, and the material sufficiently non-degenerate, to settle into quiet equilibrium burning until the 
remaining hydrogen is consumed. Multiple peaks in the luminosity evolution of ultra-luminous red 
novae are a natural consequence of this expansion- contraction cycle. 



We discern three main episodes of mass loss, spanning about 150 days, preceded by a short 
and very mild episode some 50 days earlier. During each episode the velocity rises to a maximum 
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value between 600-800 km s . The contraction of the WD at the end of each mass loss episode 
is shown by the sharp rise of the effective temperature on the one hand, and the rise in nuclear 
luminosity on the other. Expansion of the remnant shell follows, which is shown by the declining 
effective temperature (see Figure 7) and nuclear luminosity. At the end of the last mass loss episode 
we note a sharp drop of the nuclear luminosity, marking the end of the entire mass loss phase of 
the cycle. 

Again in contrast to the Paper II models, characterized by higher core temperatures and/or 
WD masses, the outburst of the novae of Table [5] occur in two stages. First, a sharp rise of the 
bolometric luminosity to about 1O 3 L lasts about one year. We predict that this year-long "pre- 
maximum" luminosity phase must occur in all ultra-luminous, red, thermonuclear-powered novae. 
Only after a year does the luminosity rise to its maximum, accompanied by mass loss. The final rise 
in luminosity - spanning a few weeks - is of only about 6 magnitudes. The decline of the bolometric 
luminosity is also very slow, typical of low mass, cold WDs. After the initial rapid expansion to 
about one solar radius noted above, the expansion proceeds at a slower rate. Significant mass loss 
commences only when the radius has increased to about W0R & , typical of nova outbursts. The 
initial expansion and rise in luminosity occur when the convective zone (starting at the base of 
the burning shell) extends all the way to the surface. The relatively slow expansion is due to the 
unusually large mass (and inertia) of the envelope as compared to the less extreme models of Papers 
I and II. 

The very high mass-loss rate (of order 10~ 3 Mq yr _1 ) associated with ejecting the massive en- 
velopes discussed here from a 0.5M & WD in just a few months would normally demand luminosities 
well in excess of 10 6 L©, if those envelopes were promptly ejected from the WD surfaces. We em- 
phasize that this is not the case. Mass loss only begins when the envelope has expanded to roughly 
lOOi?0. This bloated configuration arises because of the slow expansion phase, powered for about 
1 year by a thermonuclear luminosity of 2 x 10 5 Lq. This slow expansion phase largely unbinds the 
envelope, making subsequent mass-loss much easier to achieve with more modest luminosities. 



5. CONFRONTING MODELS WITH OBSERVATIONS OF M85-OT and M31-RV 

Soker & Tylenda (2006) have compared and contrasted classical novae, born-again AGB stars 
and stellar mergers as possible models for very luminous, red eruptive variables. They considered 
multiple characteristics in their comparisons between models and theory: increase in luminosity 
factor, multi-outburst light curve, fading as a very cool supergiant, peak luminosity, solar ejecta 
abundances, outflow velocity, association with a young B3V star, and the circumstellar non-ionized 
matter illuminated by the outburst. They concluded that the born-again AGB model fails almost 
all of these comparisons, and we concur. They also concluded that classical novae fail on six of the 
nine criteria; we disagree, and detail below how, in fact, classical novae are compatible with all of 
the observations. Fur thermore, we note that the merger model fails to explain the observed blue 



remnant for M31-RV (jShara et al,ll201ll ). whose properties are summarized below. The nova model 



naturally predicts this observation. 



5.1. Mass of nova ejecta 



The ejected masses obtained in this study are the highest that have ever been predicted by 
the classical nova TNR scenario. We have shown that 2 — 3 x 1O _3 M0 is the upper limit of m e ,- 
for classical novae. Suggestions of even larger masses - 10~ 2 M & or even KHMq - have been made 
for M31-RV, based on the assumptio n of energy equipar tition in the ejecta between kinetic energy 
and photons . All of the mod els of lYaron et al.1 (|2005l ) show that the equipartition assumption 
is incorrect (jShara et al.ll2010l ). Radiated energy dominates kinetic energy by one to four orders 
of magnitude in virtually all nova eruptions. Nova ejecta masses based on the assumption of 
equipartitio n of energy are si mply wrong, overestimating the ejected mass by at least an order of 
magnitude (jShara et al.ll2010l ). 



5.2. Nova remnant 

Another critical prediction of the nova model is that "red novae" must leave hot, blue remnant 
stars - cooling white dwarfs. We have used HST to search for such a hot, blue remnant at the site 
of M31-RV. In an accompanying paper we note the dete ction of just such an object, fading on a 



timescale of years, very close to the position of M31-RV (jShara et al.ll201ll ). That blue object was 
observed to be radiating at least 10 3 Lq when it was observed by HST seven years after M31-RV 
erupted. Twelve years after the eruption it was still radiating at 25OL , and displaying an effective 
temperature of about 40 kKelvins. Twenty years after the eruption the luminosity of the object 
was still W 2 L@, and its effective temperature no cooler than 8kKelvins. The observed cooling and 
fading behavior of this blue remnant is in perfect accord with the nova models of this paper. The 
mergeburst model predicts a much cooler (3000 Kelvins) object that fades much more slowly (on a 
timescale of many nillenia) than is observed. 

We predict that such a similar blue object must also eventually emerge from the the optically 
thinning ejecta in M85-OT if the nova model for that object is correct. 



5.3. Stellar Populations 



M31-RV erupted in the bulge of the Andromeda galaxy. iBond fe Siegell (|2006l ) examined the 
site of the eruption in Hubble Space Telescope images and reported that "there is no evidence for 
any significant young population at this location in the M31 bulge". 



M85-OT erupted near the lenticular galaxy (Hubble type SO) M85. iKulkarni et al.l (|2007l ) 



examined the site of the eruption in HST, Spitzer Space Telescope and Chandra X-ray Observatory 
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images and reported that "There is no evidence for a bright progenitor, and nor do we see tracers of 
massive star progenitors." These same authors summarized their search with "We conclude that the 
M85-OT2006-1 probably arises from a population of stars with masses of a few times the mass of 
the Sun or less" . The stellar populations around M31-RV and around M85-OT are thus consistent 
with those expected for classical novae. 



5.4. Oxygen 

A prediction of the models presented in this paper is that the ejecta in luminous red novae 
will be extremely oxygen rich, at the expense of carbon. The relatively cool peak temperatures and 
massive envelope of the low mass white dwarf we studied here allowed the CNO bi-cycle to "cook" 
about 90% of its initial C, N and O into oxygen and nitrogen. 



5.5. Carbon isotopes 

The nova models of this paper produced more Carbon-12 (relative to Carbon-13) than any we 
have ever computed. During the latter stages of ejection, C12/C13 ratios of 6 were encountered in 
the m odels. This should be compared with typical interstellar values of 60 to 90 (jMerril &: Ridgway 



19791 ). 



5.6. Luminosity increase 



An equally strong prediction of our low mass WD nova models is a pre-maximum rise to, 
and year-long plateau in luminosity at about 1000-L©. Neither M31-RV nor the M85 variable was 
observed long enough, or with sufficient sensitivity before maximum to know if it displayed this 
behavior. Multiyear synoptic surveys (su ch as the Paloma r Transient Factory (|Law et al 1 12009) 
and the Large Synoptic Survey Telescope (jAbell et al Il2009l )) will surely find more such luminous, 
red eruptive objects. We predict pre-maximum outbursts, lasting about 1 year, with luminosities 
of about IOOOLq for the future luminous red eruptive variables that are extreme classical novae. 



5.7. Peak luminosity 

A challenging observation to explain for the classical nova scenario is the high claimed luminos- 
ity for the M85 transient. At peak brightness it was about 3 times larger than the most luminous 
super-Eddington nova models we have presented in this paper. However, it is important to note 
that the nuclear luminosities of model novae during flashes - including the models calculated here 
- exceed the apparent luminosities of even M85-OT by several orders of magnitude. Modeling the 
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radiated luminosity depends critically on the expanding envelope opacity. 



The OPAL opacities (jlglesias Rogers Ill996l ) used in our simulations are only valid for tem- 
peratures greater than 6000 Kelvins. They don't include the effects of molecules that become 
important at lower temperatures. Temperatures below 6000 Kelvins have not been important in 
nova simulations until now because the much lower ejected masses of all previous models became 
optically thin with te mperatures in excess of 6000 Kelvins. Recent low temperature opacity tables 
([Ferguson et all 120051 ) include the effects of hundreds of millions of atomic and molecular lines, 
grains, negative ions, and b ound-free and free-fre e opacity. An important and remarkable result is 
summarized in Figure 11 of iFerguson et all ((2005). The Rosseland mean opacity k decreases mono- 
tonically, by five orders of magnitude, as the temperature decreases from 10,000 to 3000 Kelvins, k 
then rises (almost) monotonically by four orders of magnitude as the temperature decreases from 
3000 Kelvins to 700 Kelvins. Furthermore, At is a factor of 20 to 25 smaller at 3000 Kelvins than 
at 6000 Kelvins. This deep, well-defined opacity minimum has important consequences for nova 
envelopes - which we have encountered for the first time in the simulations presented here - that 
are massive enough to be optically thick at 3000-4000 Kelvins. 

Massive nova envelopes, like the ones in the present simulations, must "leak" large amounts of 
radiation much faster than less massive nova envelopes whose effective temperatures always remain 
in excess of 6000 Kelvins. The maximum record ed luminosities of M31-RV and M85-OT were 
8 x 1O 5 L ({Rich et al.l Il989l ) and almost 10 7 L Q (jKulkarni et al.1 12007I ). These super-Eddington 
luminosities were observed when those objects displayed the remarkably red colors of early M 
supergiants, corresponding to effective temperatures of about 4000 Kelvins. While it is beyond the 
scope of the present paper to include a realistic model atmosphere with low temperature opacities 
in the simulations of low mass WD novae, we can quantitatively estimate the highest luminosity 
that any nova can attain. For a given WD mass Mw d, the luminos ity L of a nova envelope radiating 
near the Eddington limit is given by equation 2b of IShara I (|198ll ): 



log(L/L Q ) = 4.59 + log{M wd /M G ) + log(n es /K) 

where n is the envelope opacity and n es is the opacity due to electron scattering. The key result is 
that the nova luminosity is inversely proportional to the opacity of its expanding envelope. At 3000 
to 4000 Kelvins a nova envelope's opacity (as noted above) will be 20 to 25 times smaller than the 
opacity used in the simulations described above. We therefore estimate that the highest luminosity 
attainable by a nova like the ones we have simulated (which reach 3 x 1O 5 L ) is roughly 25 times 
larger: i.e. 7.5 x 10 6 L©. This is very close to the luminosity observed for M85-OT. As we've already 
noted, there is more than sufficient energy generated (and "bottled up" ) in the model envelopes to 
power a luminosity spike to 10 7 L Q . We suggest that it is the opacity minimum noted above that 
permits the rapid leak of energy corresponding to the energy spike observed in M85-OT. 
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5.8. Fading as a very cool Supergiant and outflow velocities 

It is often noted that the very late spectral types (L or T) and cool temperatures of M31-RV 
and M85-OT are completely different from the behavior of post-outburst classical novae. This is 
because most classical novae eject shells of order 1O -5 M0 - at least a hundred times less massive 
than the ejecta of M31-RV-like objects must be to become as cool as is observed. Massive ejected 
shells (10" 3 M© ) must remain optically thick to radii 5-10 times larger than those of most classical 
novae. As noted above, their effective temperatures will cool to values considerably lower than 
those of classical nova blackbodies before they become optically thin - i.e. as low as 3000-4000 
Kelvins, corresponding to the opacity minimum noted above. As these model shells are ejected at 
velocities of a few hundred km per second (10 times lower than those of most classical novae, and 
in accord with observations), they must remain optically thick for years instead of months, just as 
is observed. 



5.9. V838 Mon association with a young B3V star 



V838 Mon lies v e ry close to (but is not co incident with) a B3V star in a young open cluster 
(jAfsar fe Bondll2007l ) . iTvlenda fc Soker I (|2006l ) emphasize that there is not enough time during the 
lifetime of a B3V star to form the white dwarf needed for a classical nova eruption, and we agree 
(barring a very rare dynamical collision event). If the erupting object V838 Mon is a companion of 
the B3V star, or if it is a co-eval member of the young open cluster then it almost certainly cannot 
be a classical nova. 

Remarkably, no such young stellar populations exist near M31-RV (and probably not M85- 
OT or M99-OT). A B star at the location of M31-RV would have been detected in archival HST 
images... it isn't there. Thus M31-RV cannot be ruled out as a classical nova. 



6. SUMMARY 

It is common to suggest that a few very luminous, red eruptive variables cannot be classical 
novae, and must therefore represent a new class of astrophysical objects. In this paper (and in an 
accompanying one) we show that models of classical novae can, in fact, reproduce the observed phe 



nome na rather well for M31-RV, M85-OT and the newly discovered red nova in M99 (jKasliwal et a] 



2010) - but not v838 Mon. The underlying white dwarfs in such nova systems must be low in mass, 
cold, and accrete at a low rate from their companions. Predictions of the nova theory are in re- 
markably good accord with observations of the "ultra-luminous red novae" M31-RV and M85-OT 
in terms of time scales and light curves. 

The key question - are "luminous red novae" really extreme classical novae or a new astrophys- 
ical phenomenon like mergebursts - will be settled when we detect more such objects well before 
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maximum light, and when we recover more post-outburst remnants after the ejecta expand and 
become optically thin. The merger of two stars must swell the resulting star, leaving behind a 
cool, red remnant. The swollen remnant's envelope must contract, becoming hotter and bluer on a 
thermal timescale (of many millenia). This predicted behavior is in sharp contrast to that expected 
from a classical nova. 

Classical novae also fade, but on a much faster timescale - years to decades. The underlying 
white dwarfs and accretion disks of old novae cool in the decades following an eruption, and thus 
become redder. The object we have detected at the site of M31-RV attained a luminosity that is 
consistent with it being a classical nova. It is observed with HST to have been very hot about a 
decade after its eruption. It is fading on a timescale of decades, and becoming much redder as it 
does so. All of these observations are in accord with the nova models of this paper, but not with 
mergeburst models. 
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Table 1. Characteristics of The Outburst - Paper II models 

Parameter Combinations Outburst Characteristics 

M WD T WD logM m acc m e j Z e j V a v Ts,max i4,mai t'A,bol 

(M ) (10 6 K) (M e yr- 1 ) (M Q ) (M ) (km s' 1 ) (10 8 K) (10 4 L e ) (days) 

0.40 10 -11 5.9E-4 7.0E-4 0.21 460 0.95 23.9 5.5E4 

0.65 10 -12 3.9E-4 6.7E-4 0.45 220 1.59 20.7 8.7E3 



Table 2. Characteristics of The Outburst - Models with M WD = O.5M 

Parameter Combinations Outburst Characteristics 

M\YD TwD M m aC c m e j Z e j V av Ts, m ax L4,,max t^.bol 

(M ) (10 6 K) (M yr" 1 ) (M Q ) (M Q ) {km s' 1 ) (10 8 K) (10 4 L ) (days) 
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1.7E-3 
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1.38 
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5E-11 


2.0E-3 


2.2E-3 


0.10 


329 


1.44 


20.8 


4.2E3 


0.50 


2 




7E-11 


1.2E-3 


1.4E-3 


0.13 


336 


1.34 


26.1 


2.7E3 


0.50 


3 




1E-10 


6.7E-4 


7.4E-4 


0.14 


104 


1.38 


17.9 
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Table 3 


. Ejecta Composition 


- Models with M WD = O.5M 






Parameter Combinations 


Mass Fractions 




Isotopic Ratios 








M WD 


Twd 


M 


12 C 


14 N 


16Q 


13 C/ 12 C 


15 N /14 N 


17 0/ 16 






0.50 


4 
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7.8E-3 


4.8E-2 
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Table 4. Characteristics of The Outburst - Very low Twd 



Parameter Combinations 






Outburst Characteristics 






M WD 


Twd 


logM 


^acc 


m e j 


Z e j 




^8, max 


-^4, max 


t'J,,bol 


(Mq) 


(10 6 K) 


(Mq yr- 1 ) 


(Mq) 


(Me) 




(km s _1 ) 


(10 8 K) 


(10 4 L Q ) 
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1.79 


-10 


9.2E-4 
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Fig. la. — Temperature profiles at the beginning of the accretion phase for 5 consecutive cycles 
for a given WD and accretion rate, as obtained in the present calculations (first model in Table 
[2]). The profiles in each panel are slightly shifted in mass with respect to one another, in order to 
enable the comparison. 
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Fig. lb. — Temperature profiles at the beginning of the accretion phase for 5 consecutive cycles 
for a given WD and accretion rate, as obtained in Paper II, for a 0.65Af Q WD of temperature 
30 x 10 6 K, which has accreted mass at a rate of 1 x 1O _11 M0 yr -1 . The profiles in each panel are 
slightly shifted in mass with respect to one another, in order to enable the comparison. 
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Fig. 2. — Nuclear luminosity history for a 0.5M Q WD of temperature 2 x 10 6 K, which has accreted 
mass at a rate of 7 x 10~ n M Q yr _1 , in the centuries before and after eruption 
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Fig. 3. — The luminosity generated in the first thermonuclear pulse of a O.5M WD of temperature 
2 x 10 6 K, which has accreted mass at a rate of 7 x 10™ 11 Mq yr _1 , a year before the optical outburst 
begins 
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Fig. 4. — The first outburst: pre-mass loss phase of an outburst cycle obtained for a 0.5M & WD 
of temperature 2 x 10 6 K, which has accreted mass at a rate of 7 x 1O -11 M yr" 1 : top left The 
temperature history at the envelope base of the first thermonuclear pulse, a year before the major 
outburst begins; top right radius of the WD envelope; bottom left envelope effective temperature; 
bottom right envelope apparent luminosity. 
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Fig. 5. — The outburst cycle obtained for a 0.5Mq WD of temperature 2 x 10 6 K, which has 
accreted mass at a rate of 7 x 1O -11 M yr _1 : top left The temperature history at the envelope 
base; top right radius of the WD envelope; bottom left six nuclear luminosity pulses generated in the 
white dwarf envelope ; bottom right envelope apparent luminosity. This peak apparent luminosity 
is artificially limited by the opacities used in the simulations, and can be 10 to 100 times larger 
when the envelope effective temperature is 3000 - 4000 Kelvins. 
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Fig. 6. — The mass loss phase of an outburst cycle obtained for a 0.5M Q WD of temperature 
2 x 10 6 K, which has accreted mass at a rate of 7 x 1O _11 M0 yr -1 : top left mass loss rate; top right 
expansion velocity; bottom left effective temperature; bottom right nuclear luminosity. 
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Fig. 7. — The outburst effective temperature obtained for a 0.5Mq WD of temperature 2 x 10 K, 
which has accreted mass at a rate of 7 x 10 yr _1 . The model is the same as that of Fig. [6] 

and time is set to zero at the same point. 



